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Abstract 

Some variants of the MSSM feature a strip in parameter space where the lightest neu- 
tralino x is identihed as the lightest supersymmetric particle (LSP), the gluino g is the 
next-to-lightest supersymmetric particle (NLSP) and is nearly degenerate with y, and the 
relic cold dark matter density is brought into the range allowed by astrophysics and cos¬ 
mology by coannihilation with the gluino NLSP. We calculate the relic density along this 
gluino coannihilation strip in the MSSM, including the effects of gluino-gluino bound states 
and initial-state Sommerfeld enhancement, and taking into account the decoupling of the 
gluino and LSP densities that occurs for large values of the squark mass niq. We hnd that 
bound-state effects can increase the maximum for which the relic cold dark matter den¬ 
sity lies within the range favoured by astrophysics and cosmology by as much as ~ 50% if 
mq/rrig = 1.1, and that the LSP may weigh up to ~ 8 TeV for a wide range of mq/rrig < 100. 
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1 Introduction 


In the absence of any signal for snpersymmetry dnring Run 1 of the LHC (^, it is natural to 
ask how and where supersymmetry may be hiding. Perhaps it is hiding in plain sight with 
a compressed spectrum that the conventional missing-energy searches at the LHC have 
been unable to resolve? Or perhaps R parity is violated, in which case supersymmetry may 
be hiding among the jets and leptons produced by Standard Model processes? Or perhaps 
R parity is conserved, but supersymmetric particles are too heavy to have been detected 
during Run 1 of the LHC? 

There are two issues with this last possibility. One is the accentuation of the problem of 
the naturalness (or hne-tuning) of the electroweak scale that low-scale supersymmetry was 
postulated to mitigate, and the other is the cosmological cold dark matter density. The cold 
dark matter may well not consist only, or even predominantly, of the lightest supersymmetric 
particle (LSP). However, even if the cold dark matter density is considered only as an upper 
limit on the relic LSP density, it imposes an upper bound on the LSP mass that depends on 
the specihc LSP candidate under consideration. 

If R parity is conserved and the LSP is present in the Universe today as a relic from the 
Big Bang, it is expected to be electromagnetically neutral and have only weak interactions. 
In the minimal supersymmetric extension of the Standard Model (MSSM), the most plausible 
candidates are the lightest neutralino y and the gravitino [^. Here we consider the neutralino 
case, and the cosmological upper bound on its mass. 

The relic LSP density depends not only on the LSP mass, but also on the rates at which 
it annihilated with itself and coannihilated with other sparticles in the early Universe [^. 
Other things being equal, the largest LSP mass is allowed when such coannihilation rates 
are maximised, which happens when the LSP is (nearly) degenerate with other particles. If 
there is only one such coannihilating sparticle species, the coannihilation rate will in general 
be maximised for a coloured sparticle. There have been analyses in the literature of the cases 
where the coannihilating particle is a squark, specihcally the lighter stop squark [5 11 , and 
also the case of the gluino ^ 12 -18 . In general, one would expect that the heaviest LSP 
will be allowed when it coannihilates with the particle with the largest colour charge, namely 
the gluino. 

We study here the question how heavy the neutralino LSP y could be, if it is nearly 
degenerate with, and coannihilates with, the gluino g. This is of relevance to assessing, 
for example, what centre-of-mass energy would be needed for a proton-proton collider to 
be ‘guaranteed’ to detect R-conserving supersymmetry. There can of course be no cast- 
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iron guarantee, even within the MSSM. For example, even in the gluino coannihilation case 
studied here the neutralino LSP mass limit depends on the squark masses, and the LSP mass 
limit could be substantially modified if the squarks were degenerate with the neutralino LSP 
and the gluino. However, a complete analysis of this case lies beyond the scope of this paper. 

As already mentioned, there have been several previous analyses of neutralino-gluino 
coannihilation 10,12 -18 , and the main new elements here are in our discussions of the 


effects of gluino-gluino bound states and of the issue whether coannihilations can be main¬ 
tained in the presence of a large squark to gluino mass ratio. Here, we will restrict our 
attention to the coannihilation processes and leave their application to more complete mod¬ 
els (with for example, radiative electroweak symmetry breaking) for future work 19 . As 
we discuss in detail, bound states can remove from the primordial plasma gluino pairs that 
may subsequently annihilate into Standard Model particles, before they can decay into the 
LSP as is usually assumed in discussions of coannihilation. We present numerical estimates 
of the bound-state production rate, and find that, for fixed sparticle masses, the relic dark 
matter density is substantially reduced compared with the cases where bound-state forma¬ 
tion is neglected. Conversely, the cosmological relic density may lie within the cosmological 
range for substantially larger LSP masses than would have been estimated in the absence of 
bound-state effects: this effect is ~ 50% for mg/rrig = 1.1, falling to ~ 20% for mq/rrig ~ 10 
to 50. Another effect we discuss is that, if mg/rrig > 100, the densities of neutralinos and 
gluinos decouple and coannihilation effects freeze-out early, leaving a significantly higher relic 
density, thereby reducing the possible LSP mass. There is also a reduction in the possible 
LSP mass for small rriq/mg —)■ 1, due to cancellations between s-, t- and n-channel diagrams 
that tend to reduce annihilation rates. 

Taking these effects into account, we find a maximum LSP mass ~ 8 TeV if it is the 
Bino, which may be attained for 10 < triq/mg < 100. If the LSP is the neutral Wino, the 
upper limit is reduced to ~ 7 TeV, and for a neutral Higgsino the upper limit it becomes 
~ 6 TeV. 

The layout of this paper is as follows. In section 2, we review the Sommerfeld en¬ 
hancement in the relevant gluino-gluino annihilation processes. In Section 3, we discuss the 
formation of gluino-gluino bound states, considering also dissociation processes in the early 
Universe. In Section 4, we consider the rates for conversion between gluinos and neutrali¬ 
nos. In Section 5, we present and discuss the coupled Boltzmann equations for neutralinos 
X, gluinos and gluino-gluino bound states. Section 6 contains some numerical results for 
the gluino coannihilation strip and a discussion of its endpoint. Section 7 summarises our 
conclusions and discusses their significance for future colliders. Finally, Appendices present 
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technical aspects of the computation of the 2 —)■ 2 cross sections needed for solving the 
Boltzmann equations. 


2 Sommerfeld Enhancement 


Before discussing the formation and effects of gluino-gluino bound states, we hrst discuss 
briefly Sommerfeld effects in gluino-gluino annihilation, which may enhance annihilation 
rates at low velocities, and are particularly relevant in the case of the strongly-interacting 
gluino. As a general rule, initial-state interactions modify s-wave cross-sections by factors 


20,21 


F{s) ^ 
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a 


( 1 ) 


where a is the coefficient of a Coulomb-like potential whose sign convention is such that 
the attractive case has a < 0, and f3 is the velocity of one of the annihilating particles in 
the centre-of-mass frame of the collision. In the cases of strongly-interacting particles, the 
Coulomb-like potential has the form (2^ 


CXa 




( 2 ) 


where as is the strong coupling strength, C/ is the quadratic Casimir coefficient of a specific 
final-state colour representation, and Q and C[ are the quadratic Casimir coefficients of the 
annihilating coloured particles. In the case of octet annihilating particles such as gluinos, 
Ci = C[ = Cg = 3. The relevant hnal states are in singlet, octet, or 27^ representations, for 
which Cf = 0, C/ = 3, or Cf = 8. 

As discussed in [^, Sommerfeld effects such as these have been implemented in the SSARD 
code 23 for calculating the relic dark matter density. In the coannihilation region of interest, 
this code uses a non-relativistic expansion for annihilation cross-sections: 


{av) = a + hx ^ -f ... 


( 3 ) 


where (...) denotes an average over the thermal distributions of the annihilating particles, the 
coefficient a and h represent the contributions of the s- and p-wave cross-sections, x = m/T, 
and the dots represent terms of higher order in 1/x. A Sommerfeld enhancement occurs when 
a < 0 in Q, modifying the leading term in (|^ so that it acquires a singularity oc ^/x. In 
this paper we have included these enhancements in the gg —)■ gg and gg —qq cross sections. 
The procedure for obtaining a thermally averaged cross section is given in Appendix A. The 
expressions for the matrix elements for the coannihilation processes are given in detail in 
Appendix B. 


3 







3 Gluino-Gluino Bound-State Formation 


Gluino-neutralino coannihilations may increase the effective annihilation cross section and 
thereby lower the final nentralino relic abnndance. The Sommerfeld enhancement discnssed 
above farther increases the cross section in specific channels and again lowers the abnndance 
of nentralinos allowing for larger masses at the tip of the coannihilation strip defined by 


Am = 0 where Am is the glnino-nentralino mass difference 17 . Glnino-glnino bonnd states 
can farther serve to remove glninos from the thermal bath and thereby lower the relic density 
by a factor that is non-negligible relative to the Sommerfeld enhancement, and mnch larger 
than the nncertainty in the cosmological cold dark matter density. 

The dominant process for the formation and dissociation of glnino-glnino bonnd states 
R in the thermal plasma is g + g R + g. These processes become important when the 
plasma temperatnre falls low enongh for typical thermal energies to become comparable to 
the binding energy of the R state, namely T < Eb = 2mg — In principle, one may form 
colonr-octet states as well as singlets, bnt the latter are expected to be more deeply bonnd 
with larger wave fnnctions at the origin. Here we focns on the prodnction of the lightest 
colonr-singlet state, 1^, with orbital angnlar momentnm L = 0 and spin angnlar momentnm 
S' = 0, which is expected to be the most copionsly prodnced. Since we are considering 
glninos weighing several TeV, we expect the leading order of QGD pertnrbation theory to be 
a nsefnl approximation, and assnme the Gonlomb potential V{r) = —for the state, 
with binding energy Eb — {3as/2)^mg. The normalised spatial part of the wave fnnction for 
this lo bonnd state is 


^(r) = 


( 4 ) 


where a = 2/{3asmg) is the Bohr radins. The bonnd state decays predominantly to a 
pair of glnons and the leading order decay rate is 

243 . 


~ ^ ot^mg. 


( 5 ) 


3.1 Dissociation 


In order to calculate bound-state formation and dissociation via the dominant processes 
g°‘ + g^ R + g^, we first calculate the bound-state dissociation cross section, adis following 
Section 56 of 24 , where the photoelectric effect for an atom is calculated. The central part 


of the calculation is the evaluation of the transition amplitude given in Eq. (56.2) of 24 

f V ■ ? - 

Mh = 




( 6 ) 


4 






where 0/ is the wave function of the free g°' pair and 0* = 0bs(r), and and k are the 
polarisation and momentum vectors of the gluon, respectively. 

We use the dipole approximation, ~ 1, which is justihed because the bound-state 
wave function 0bs(r) is exponentially suppressed for r > a, and because the gluon momentum 
\k\ = u, where its energy u satishes the conservation condition 

a;2 |pl2 

ZyZTflg ^B) ^^9 

where |p| is the momentum of one of the annihilating gluinos. (Note that |p| is the same as 
the relative momentum, {mg/2)vreu and the second term on the LHS of ([^ can be neglected.) 
We hnd uja ~ Esa = <C 1 for Vrei = 0 and = 0.1, and more generally ua < 1 when 
Vrei < 0.6, so that the dipole approximation should be sufficient for our purposes. 

The dipole approximation imposes a selection rule on 0/, which needs to be in an 
L = 1 state. Further, charge conjugation (C-parity) conservation requires that C{g°‘g^) = 
C{R)C{g'^), where the ground state with L = 0 and S' = 0 has = 0“+. The C- 
parity of the colour anti-symmetric 8^ state is the same as that of the gluon 25 , while the 
C-parity of colour-symmetric 8^ state is opposite of that of the gluon, for all color indices. 
Therefore, the only possible state for (j)f is 8^, with L = 1 and S' = 0. (Note also that parity 
is conserved in this case, because P{(j)f) = 1 and the gluon has P = —1.) 

The normalised spatial part of the wave function for the free pair g°'g^ is 


- OO 


( 8 ) 


Only the L = 1 term survives, due to the selection rule from the dipole approximation. Since 
we wish to calculate \A4fi\‘^, we may discard the phase shift factor is real) and the 

factor Therefore we write 
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where A4 ji is calculated following Section 56 of . 


Since 0/ is the wave function for an 8^ state, the Coulomb potential is V/(r) = — |as/r. 


whereas 0* is a wave function for the Coulomb potential V (r) = —3as/r, the result is different 

^ -|-^^) g-4?arctan(l/K^) ^ 


from Eq. (56.12) of 24 , namely 
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where ^ = ^ag/vrei and k = 2. This equation is averaged over the gluon polarizations and 
would reduce to Eq. (56.12) if k = 1. 

The total wave functions for the free pair and the bound state R are products of 
the spin, colour and spatial parts of the wave functions. In view of the Majorana nature of 
the gluinos, the total wave functions should be anti-symmetric. Concerning the spin part of 
the wave function, since both the bound state and the free gluino pair are in an S' = 0 state, 
the spin wave functions are both 

(12) 


and the spin parts of the wave functions do not introduce any extra factor in adis- As for the 
colour part of the wave function, according to 26,27 (pi,color = Sde/ \/8, and (pf,color = fhacj \/3 
(the latter because fabcfabd = 3(5cd)- 

The (~^ (m-^ 2 ) )^^^ ^ factor in the transition amplitude (6) is calculated from the gluino- 
gluino-gluon interaction Lagrangian 




(13) 


which can be compared with the corresponding QED interaction Lagrangian 


C = -eQfA.h^f . (14) 

We simply replace the electric charge factor Qj in the transition amplitude (|^ by the colour 
factor fabc, since the factor 1/2 in (13) is compensated by a factor of 2 due to the Majorana 
nature of the gluino. Putting the above colour factors together, we obtain 


1 fhad n |2 _ 1 ^ |2 _ o 


(15) 


Vs ^3 

Note that all color indices are summed over. 

Concerning the spatial part of the wave function, we need to take into account the fact 
that both the initial and hnal states contain two identical particles. In the case of the bound 
state, they are in the symmetric L = 0 state, and the wave function needs to be symmetrised 


as in Eq. (2.14) of 26 


[<i>bs(ry 9, 4b} + 4bsir, r - $,4 + Tr}]^ ^2 X (;ra^) 


(16) 


On the other hand, the hnal free pair is in the antisymmetric L = 1 state, and the wave 
function needs to be antisymmetrised: 


,p ■ r ^ 


— [(pf{r,e,(p)-(pf{r,TT-e,(p + 'K)\ = \/2 X 


(17) 
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The coefficients in these two equations introduce an extra factor of |a/ 2\/2P = 4 into the 
modulus-squared of the spatial wave function factors. Finally, recall that we have averaged 
over the polarisations of the gluon, but we also need to average over its colour. This gives a 
factor of 1/8. Therefore, the final dissociation cross section is 

adis = 3x4x^xix al ,, (18) 

where the hnal factor of 1/2 is to avoid double counting of gluinos in the hnal-state phase- 
space integration. 


3.2 Formation 


We come hnally to the bound-state formation cross section, absf, which is related to adis 
through the Milne relation: 

^nfnfabsfVrei (^1 + ^uj/t _ x ) f\'*^rei)dVrei = rf^adisdn^"^ , (19) 

where the | on the LHS of the above equation is introduced to avoid double-counting the 
number of bound-state formation reactions, and the factor ,,T ^ comes from the enhance- 
ment of bound-state formation due to the stimulated gluon emission in the thermal back¬ 
ground (similar to the stimulated recombination in e~p O TTy). Using 


dnl^ 

f i^rel) 


An uj'^doj 


9g 


(27r)^ 


1 ’ 




2nT 


2 ,g-(m9/2)Kei 
rel^ 


/2T) 


( 20 ) 


and ([^, we hnd 


where 




gpQg _ 1 X (2 X 8) ^ ^ 
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For comparison, the Sommerfeld enhanced s-wave cross section for —)■ g'^g’^ is given in 

Eqs. (2.13) and (2.25) of 


. . 27r(2e) 

^ann\S^annVrel) I ^ 


^ 1 27re ^1 27r(-|0 ^ [Anas? 27 


6 1 _ e-2^(2«) 3 1 _ e-2^C 2 i _ J ml 5127r ’ 


(23) 


7 













where ^ = ^ag/vrei- Therefore, in the Vrei —t 0 limit we hnd 

^ann\S^an7i^rel) 

Therefore, we see that the inclusion of gg bound states is a non-negligible component in 
determining the hnal neutralino relic density. 

4 Conversion Rates 

For coannihilation to be effective, the coannihilating species (in this case neutralinos and 
gluinos) must be in thermal contact. That is, the rates for interconverting the LSP and 
NLSP must be faster than the Hubble rate. In both the familiar cases of stop and stau 
coannihilation, connectivity of the two species can be taken for granted, as the conversion 
rates are mediated by light Standard Model particles and are always fast. This implies that 
the ratio of densities {riNLSp/nLSp) is approximately equal to the equilibrium ratio and allows 
for a simplihcation in the Boltzmann equations. However, the interconversion of neutralinos 
and gluinos must proceed via squarks, leading to a suppression if the squarks are heavy. 
The relevance of the coannihilation process relies on fast conversion rates, and requires the 
ratio of squark masses to the gluino mass to be less than approximately 100 as we show 
below. For larger squark masses, the gluino and neutralino abundances evolve separately, 
and coannihilation effects are essentially shut off independent of the mass difference. 

The interconversion processes we consider are X 5 ' -H- yg O gg, and the gluino decays 
and the inverse decays g <H- XQQ- When the neutralino is a Wino or Higgsino, the processes 
involving a chargino, x^d <H- gu, x'^u -H- gd and g -H- x~^du, as well as the corresponding 
processes for x~, are also included. We note that g stands here for all six quark flavors, and 
the u, d stand for all the three generations of up-type and down-type quarks. Also, when y 
is a Higgsino, the two lightest mass-degenerate neutralino components, Hi^ 2 , are both taken 
into account. For each relevant process, we first calculate the transition matrix element |Tp. 

We calculate the gluino decay rates for g —)■ ygg, g —)■ x^du and its charge-conjugated 
process]^ The squared transition matrix elements |Tp are identical to the corresponding 

^Since we treat quarks as massless, three-body decays will always occur as long as rrig > m^. Due to 
logarithmic corrections, the gluino two-body decay into a neutralino and a gluon becomes important when 
the squark-to-gluino mass ratio is very large 100) [28| . However, as we show below in Section 6, gluino 
coannihilation becomes irrelevant for such a large squark-to-gluino mass ratio. On the other hand, when 
the squark-to-gluino mass ratio is ^ 100, the gluino-neutralino conversion rate dominates over the gluino 
two-body decay rate when the g and y are so close in mass that all the gluino three-body decays would be 
kinematically forbidden, if quarks are not treated as massless. 




ones for the coannihilation processes given in Appendix B, except that the expressions in 
Appendix B should be multiplied by a factor of 2, because the statistical factor for the initial 
spin averaging is i X H I for the coannihilation processes, whereas it is | for the gluino 
decay processes. We note also that the dehnitions of the Mandelstam variables should also 
be changed correspondingly as follows: for the coannihilations, s = {pi +^ 2 )^, t = {pi — ps)^ 
and u = {pi — PiY, whereas for the gluino decays, s = {pi — P 2 Y, while t and u do not 
change. 

The gluino decay rates are then obtained by performing the standard 3-body phase space 
integration. The inverse-decay processes do not have to be calculated separately, because 
they are taken into account automatically by the Boltzmann equations given in the next 
section. 

To calculate the conversion rates for XQ gq, x^d —)■ gu, x^u —)• gd and their charge- 
conjugated processes, we hrst calculate the cross sections. Again, the squared transition 
matrix elements |Tp are identical to the corresponding ones for the coannihilation processes 
given in Appendix B, except that the expressions in Appendix B should be multiplied by 
a factor of |, because the factor for the initial color averaging is | for the coannihilations, 
whereas it is | for the conversions. Also, compared to the coannihilation processes, the 
Mandelstam variables for the conversion processes are re-de£ned as s = {pi — P 2 Y, t = 
(pi ± psY and M = (pi =F PaY^ where the upper signs in the dehnition of t and u apply if qb 
or ds is brought into the initial state, while the lower signs apply if qa or ua is pulled over 
to the initial state. 

For each of the quark flavors, the thermally-averaged conversion rate is obtained by 
integrating acVq over the Fermi-Dirac distribution of the quark in the initial state. 


(rc) = 


f + OO 


apvqdnq = 


CTcVq 


3 ■ 2 ■ 47r \Pq\‘^d\pq\ 
(27r)^ -|- 1 ’ 


(25) 


where dc is the conversion cross section for any of the relevant processes discussed above. In 
the initial neutralino or chargino rest frame, cXc is a function of the incoming quark energy 
Eq. In this reference frame, f or u is the squared center-of-mass energy, and is given by 
-|- rUg -|- 2my.Eq, where the lower limit of Eq is Eq^.^ = [{rrig + mq/Y — — mq]/(2m^), 

where q' represents the quark in the hnal state. Here Vq is the velocity of the incoming quark, 
and it is related to the energy and 3-momentum of the quark by Vq = \pq\/Eq. The factors 
3 and 2 in (25) count the quark color and spin degrees of freedom, respectively. Again, the 
inverse conversion rates are taken into account automatically by the Boltzmann equations. 
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5 Boltzmann Equations 


We are now in a position to put all of the components discussed above into a rate equation 
(or set of equations) in order to solve for the relic density. To do so, we begin by considering 
three separate density components: neutralinos, gluinos and bound states. 

To set up a coupled set of Boltzmann equations, it is convenient to rescale the number 
densities of neutralinos, gluinos and bound states by the entropy density. 


Ti 

V — hr V — 11^ V- — 

5 5 ^ • 
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These are governed by the following set of coupled Boltzmann equations 


dx 


xs 


H{m^) 
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^ {r=> ( (u - yy^j + {r)i (u 




g 


dx 


xs 


H{m 


XJ " 


(UU - Ypvp) - (nu - Y‘%’) 


g 


Y, 


+ 2. (rc) y (>^x - ) - (r),y (- yp 


-Mt 


V-V — 

^g^g ^g ^g yeq 
R 




g 


dY, 


R 


XS 


dx 


H(m^ 


(Os-lUi 


yfQ 

R 




yp? 

R 


where 
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X = s = ^g^sT'\ H{m^) = H{T)x^ = 
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(26) 


(27) 


(28) 


(29) 


(30) 


T 45 

and g^s and g^ are the total numbers of effectively massless degrees of freedom associated with 
the entropy density and the energy density, respectively, is the relative velocity times 

the total cross section for the channels for yy annihilation into Standard Model particles, and 
(cTu)^- and (o'v)-- are to be understood similarly, (Tc) and (T)- are the total conversion 
rate and gluino decay rate discussed in the previous section, and all possible quark and 


set of coupled Boltzmann equations for the photino and a gluino it-hadron was studied in [^. In a 
different context, Boltzmann equations involving a bound state can be found, for example, in 


30 
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anti-quark channels for the x are summed over, (F)^ is the decay rate of the i?, and 
is the bound-state formation cross section times the relative velocity of the two incoming 
gluinos, taking into account the l/(e‘^/^ — 1) enhancement factor as discussed in Section 3.2. 
Finally, {crv )has the same effect as (r)^/s, namely, it converts the bound states to 
gluons without altering the density of free gluinos. All the quantities bracketed by (...) are 
thermally averaged, and the superscript ‘eq’ denotes equilibrium yields. 


Eq. (29) can be written in a more intuitive form: 


(r)^+ 


din a; 


H{T) 



Uav) 


n- I S 
bsf'h I V. 


H{T) 



31) 

One can check that the LHS of Eq. (31) is of order -10, whereas each of the terms on the RHS 


of Eq. (31) are of order where Mp = . Hence, to a good approximation, we 

can set the two terms on the RHS equal to each other and solve for ( ^ 

V ^ R 


Xr 

yeq 

R 


X)r + + (F)^.^ 

X)r + XX-gR^gg^g + X) dis 


where 

Therefore, we hnd that 

d(Y, + y,) 
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(32) 


(33) 


xs 


dx 


H(m,) L 


- E w>«(ry;-yyyr) 

( r>H + (^<'>8 




-{av) 


gR^gg^d 


bsf 


{r)s + + <r)8i. 


(ysYi-Yiy) ,( 34 ) 


Moreover, we note that Xv)gp^ggRg is much smaller than (F)^ for x > 10, due to the 
fact that rig decreases with the decrease of temperature while (F)^ is nearly temperature 
independent. Since the process gR —)■ gg is related to the bound-state formation process by 
crossing, should be related to {crv)^^^ by a coefficient not too much different from 

order 1. 

If at least one of the Xq Xc) and (F)- is sufficiently larger than H{T) throughout the 
period during which {Y^ + Yg) changes substantially, which is the case when the squark mass 
appearing in the denominators of the matrix elements for these processes is not too large. 


Eq. (34) can be solved by using the very good approximation YgjY^ Y^ jY^. In this case. 
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Eq. (|34|) can be recast in the familiar form snitable for coannihilation calcnlations, and we 

X dg^ 


can write 


dY 

dx 


xs 


1 - 




(r" - (K,)") , 


(35) 


V ‘^9*s dx 

where we have inclnded the term which takes into acconnt the evolntion of g^s with 

temperatnre. As we will see, this approximation is valid so long as mg/rrig < 20. 


In Eq. (35), Y = n/s, where n is interpreted as the total nnmber density, 

n = rii = n-^ + rig . 

i 

and Yeq = rieq/s, where rieq is the total eqnilibrium nnmber density. 


"^eq — 


n, 


i 

The effective annihilation cross section is 

(^eff^rel) = ^ ^ 


eq ,2 


= + n 


eq 


^eq,2^eq,ji 
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(T^ji'Urel) • 


(36) 


(37) 


(38) 




eq 


As one can see from Eq. (34), the expression for {(Tv)gg is the ‘standard’ term in the hrst 


line of (34) combined with the second line involving the bonnd states. We re-emphasise that 


this simplihcation reqnires a fast interconversion rate as discnssed in the previous Section, 
so that we can set {Yg/Y-^ = {Y^/Y^^), which is true only when mg/rrig < 20. For larger 
squark masses, we use the coupled set of Boltzmann equations to solve for the relic density. 

When the LSP is a Wino or a Higgsino, we can still use all the above equations to solve 
for the relic density. All we need to do is re-de£ne the following quantities to include the 
contributions from each of the y components, y*, neutral or charged, as 



III 

M 


Xi 

"'x 

III 

M 

s 


Xi 

^^)xx 

= (^^)x^x/xYx, 


Xi,Xj 




Xi 

(Fc) 



Xi 




(39) 


Xi 
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where q' is the same as q when Xi is a neutralino, and they are different when Xi is a chargino, 
and the q and q' indicate all the possible quark and anti-quark channels for the conversion 


rates and gluino decay rates. In Eq. (39), = ^xil'^x^ where the latter ‘=’ is 

guaranteed by the fast conversion and/or decay rates among the different y/s. For later 
discussion, it is useful to dehne an effective number of degrees of freedom for y: 


3x.sf = ^9xi (1 + ^xif^ ^'^v{-^x^^xJT) , (40) 

Xi 

where A^- = — 1), and we assume Xi is the lightest component (i.e., the LSP) 

among the y/s. We can then write explicitly as 

W. = (1 + ^xJ^'^^exp(-A^,mxi/T). (41) 

9Xeff 

In the limit that all the y components have the same mass, 5'xe// = 2, 6 and 8 for Bino, 
Wino and Higgsino, respectively. 


6 Numerical Results 

We now present some numerical results obtained using the above formalism. Our results in 
this section are based on simplihed supersymmetric spectra dehned at the weak scale. We 
assume degenerate squark masses, rrig and for the most part, our results do not depend on 
supersymmetric parameters such as /r, Aq, and tan/^ We assume that the neutralino is 
a pure state of either a Bino, Wino, or Higgsino. Thus our free parameters are simply the 
neutralino mass, m^, the gluino mass, rrig and the squark masses, rrig. In future work, we 
apply these results to more realistic CMSSM-like models (without gaugino mass universality) 
and pure gravity mediation models with vector-like multiplets [I^ . 

We begin with the case in which the lightest neutralino y is the Bino. 

6.1 Bino LSP 

Fig. ID compares a naive calculation using a single Boltzmann equation for the total relic 
abundance (left panel) with a treatment of the three coupled Boltzmann equations for the 
gluino, Bino and gluino-gluino bound state abundances (right panel). These results are 
for the representative case = 7TeV, Am = rrig — = 40 GeV and mq/rrig = 10. 

The dashed line in the left panel shows the total relic abundance as would be given by the 

^The exception is the case of Higgsino-gluino coannihilations for which the vertices do depend on tan f3. 
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m^/T 


Figure 1: Left panel: The evolution of the total supersymmetric particle abundance n/s as a 
function of m^/T for the representative case = 7 TeV, Am = mg — = 40 GeV, and 

mg/mg = 10 using the single Boltzmann eguation (3f). The dashed line exhibits the naive 
thermal eguilibrium abundance, and the solid line shows the numerical solution of eguation 
(34) for the sum of the neutralino and gluino densities, exhibiting the familiar freeze-out 
when m/T ~ 30. Right panel: the separate evolutions of the abundances (solid lines) using 
the full set of eguations given by Egs. (27 - 29) for the gluino (blue lines), Bino (red lines), 
gluino-gluino bound-state (green lines) and sum of the Bino and gluino densities (black lines) 
as functions of m^/T. Details of the evolutions of the abundances are shown in the inset. 
The dashed lines again exhibit the naive thermal eguilibrium abundances. 


Boltzmann distribution if thermal equilibrium were maintained, and we see a clear departure 
for m/T > 30, as expected from a freeze-out calculation. In the right panel of Fig. the blue 
line shows the evolution of the gluino abundance, the red line that of the Bino abundance, 
the green line that of the gluino-gluino bound states, and the black line that of the sum 
of the gluino and Bino densities. The inset shows details of the evolutions of the gluon, 
Bino and bound-state abundances. The dashed lines again show the corresponding naive 
thermal equilibrium abundances. With the stated choices of m^, m^ and m^, the relic dark 
matter density using the single Boltzmann equation is = 0.120, and the full set of 

three Boltzmann equations yields = 0.119. Correspondingly, the black lines in the 

right panel of Fig. [^are indistinguishable from the lines in the left panel. 

Fig. I shows the effect of varying mq/mg for the same representative values m^^ = 7 TeV 
and Am = mg — m^ = 40 GeV. The left panel is for m^/mg = 1.1in which case we hnd 
that the relic cold dark matter density is higher than previously: = 0.21. This is due 

to the fact that at a low squark to gluino mass ratio, there is a cancellation among the t 

^We do not show results for smaller values of mq/rrig, since then squark-Bino coannihilations should also 
be taken into account. 
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Figure 2: ^45 in the right panel of Fig. but for the choices mg/irig = 1.1 (left), 120 (right). 


and u channel annihilations with the s channel leading to a smaller gluino annihilation cross 
section and hence a larger relic density. The results also change, even more signihcantly, for 
large values of mq/irig, as shown in the right panel of Fig. where mq/irig = 120. In this 
case, we hnd a much larger value of ~ ^.0. In both panels, the insets show details of the 
evolutions of the gluino, Bino and bound-state abundances. At still larger values of mg/mg 
the relic density grows very sharply: for example, for mg/irtg = 150, we hnd = 780, and 
for mq/mg = 200 we hnd ~ ^ assuming the same gluino-Bino mass diherence 

of 40 GeV. These large numbers rehect the failure of gluino-neutralino conversion to keep 
pace with the Hubble expansion for large mq/mg. 

In order to summarize the ehects of both the cancellations in the annihilation cross 
section at low mq/mg and the decoupling of the gluino coannihilations at high mq/mg, we 
show in Fig. the relic neutralino density as a function of m^/mg for our nominal value 
of m^ = 7TeV, and Am = m^ — m^ = 0,40, and 120 GeV (black, red, and blue lines, 
respectively). We see clearly the rise in Q^h^ at small m^/mg as well as the very rapid rise in 
at high m^/mg > 100. In between there is a plateau with lower Q^h‘^, as exemplihed by 
the case m^/mg = 10 shown in Fig. In general, there is a shallow minimum in around 
mq/mg ~ 50 whose location depends on Am. The horizontal band indicates the 3 -ct range 


for the Planck determination of the cold dark matter density of = 0.1193 ± 0.0014 31 


The panels of Fig. [^display bands in the (m^. Am) plane where 0.1151 < < 0.1235 

(3 a below and above the current central value [^) for a selection of values of mq/mq, as 
calculated in various dynamical approximations. The red bands were calculated dropping 
both the Sommerfeld enhancement factor and the effect of gluino bound-state formation. As 
was already noted in [9,17 the Sommerfeld enhancement causes a signihcant suppression of 
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Figure 3: The relic cold dark matter density as a function of mq/mg for = 7 TeV 
and the choics Am = mq — m^ = 0,40, and 120 GeV (from bottom to top, black, red, and blue 
lines, respectively). The rise at small mq/mq is due to the cancellations between the s-, t- and 
u-channel diagrams for gluino pair annihilation, and the rise at large mq/mq is due to the 
decoupling of the gluino and neutralino densities. The horizontal band indicates the 3-a range 
for the Planck determination of the cold dark matter density of klh? = 0.1193 ± 0.0014 (311. 


for fixed values of the model parameters, so long as the gluino-Bino conversion rates are 
large enough that the gluino coannihilation is effective. Correspondingly, the orange klydP 
bands, calculated including the Sommerfeld factor, appear at larger values of Am and extend 
to larger values of m^. The effect of including bound-state effects is to suppress further the 
value of Vtyjp for hxed model parameters, so that the corresponding black bands in 

Fig. [^extend to even larger values of Am and m^. We also show in Fig. (coloured purple) 
the bands that would be found if the bound-state formation rate were a factor 2 larger than 
our calculations, as might arise from higher-order QCD or other effects 


®We evaluate the as appearing in the Sommerfeld enhancement factor at a scale Orrig that is typical 
of the momentum transfer of the soft-gluon exchanges responsible for the Sommerfeld effect 32 , and take 


P = 0.3, which is comparable to the thermal velocities of the gluinos at the freeze-out temperature. Because 
the Sommerfeld enhancement is a precursor to the formation of bound states, for simplicity we take the 
same in evaluating the bound-state effects. The full QCD potential and the thermal mass of the gluon 
were considered in computing the Sommerfeld enhancement in 17 , and they should be relevant also to the 


calculation of the bound-state effects. We estimate that these effects could result in a shift of our orange 
bands inward by ^ 10%, based on a comparison of the light green and dark green bands in the lower right 
panel of Fig. 2 in 17 , and we expect a similar shift for our black bands. However, these two additional 


effects could be compensated by contributions from excited bound states that we do not include. We note 
that it was estimated in 33 in the context of a hidden-sector stau that excited bound states can at most 


introduce a factor of ~ 1.6 in the thermally-averaged bound-state formation cross section, compared to that 
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Figure 4: The Am = mg — m^) planes for a Bino LSP, exhibiting bands where 
0.1151 < < 0.1235 (3 a below and above the current central value), for different values 

of mq/mg = 1.1 (upper left), 10 (upper right), 50 (lower left) and 120 (lower right). These 
results are calculated without the Sommerfeld enhancement factor and gluino bound-state for¬ 
mation (red bands), with the Sommerfeld enhancement factor but without gluino bound-state 
formation (orange bands), with both the Sommerfeld enhancement factor and gluino bound- 
state formation (black bands), and allowing for the possibility that the bound-state formation 
rate is a factor 2 larger than our calculations (purple bands). 
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The upper left panel of Fig. [^is for the case mq/rrig = 1.1, where the t and u channels 
partially cancel the s-channel contributions to the gluino annihilation cross section. Here 
we see that the black band calculated including both the Sommerfeld enhancement and 
gluino bound-state effects extends to ~ 6.2 to 6.4 TeV. In this case, the numerical effects 
of the Sommerfeld enhancement are similar to those of gluino bound-state formation, and 
both effects are considerably larger than the current observational uncertainties in the dark 
matter density represented by the breadths of the bands. The purple band, which includes 
an allowance of a factor 2 uncertainty in the bound-state effects, as might arise from higher- 
order QCD, excited bound states, etc., extends to larger ~ 7.2 to 7.5 TeV. In the case 
mg/Trig = 10 (upper right panel of Fig. |^, the effect of bound-state formation is somewhat 
smaller than the Sommerfeld effect, and the black (purple) band extends to ~ 8 (9) TeV. 
These trends are also seen in the case mg/rrig = 50 (lower left panel of Fig. |^, where the 
black and purple bands also extend to ~ 8 (9) TeV. On the other hand, the results for 
mg/Trig = 120 (lower right panel of Fig. are quite different. The Sommerfeld effect is much 
larger than the bound-state effect though the latter is still slightly larger than the widths of 
the coloured bands corresponding to the S-a ranges for the cold dark matter density. Also, 
the allowed range of the LSP mass is greatly reduced, extending only to ~ 6.1 TeV (~ 6.3 TeV 
allowing for a factor 2 uncertainty in the bound-state effects). 

Fig. displays these effects differently, exhibiting the positions of the endpoints {rrig = 
m^) of the gluino coannihilation strips as functions of the assumed value of again from 

calculations with neither the Sommerfeld enhancement nor gluino bound states (red lines), 
with the Sommerfeld enhancement but without the bound states (orange lines), with both 
effects included (black lines), and allowing for a factor 2 uncertainty in the bound-state effects 
(purple lines). The horizontal green bands again show the 3-cr band 0.1151 < < 0.1235. 

We recall that lower values of would be relevant if the LSP provides only part of the 
dark matter density, e.g., if there is also a signihcant axion component, and parameter choices 
yielding higher values of in conventional Big Bang cosmology (as assumed here) could 
be relevant in models with non-standard cosmological evolution As was also seen in Fig. 
the smallest value of mg/rrig = 1.1 (upper left) leads to smaller values of for any hxed 
value of as compared to the mg/rrig = 10 case (upper right). The choice mg/rrig = 50 

obtained by considering only the ground state, but the color charge in the gluino case may change this value. 
We therefore plot purple bands with an uncertainty of a factor 2 to allow for these uncertainties. 

®It is worth recalling that for rrig ~ 5 TeV the bound-state binding energy is 0(50) GeV and the freeze- 
out temperature in conventional Big Bang cosmology is hundreds of GeV, so the assumption made in this 
paper of standard cosmological evolution is rather different from that made in more conventional thermal 
dark matter scenarios where the freeze-out temperature may be in the GeV range. 
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Figure 5: The locations of the endpoints of the gluino coannihilation strips for different values 
of Vt^h^, using the same colour conventions as in Fig. As in that Figure, the calculations 
assume the different values mq/mg = 1.1 (upper left), 10 (upper right), 50 (lower left) and 
120 (lower right). The horizontal green bands show the 3-a band 0.1151 < Q^h"^ < 0.1235. 

(lower left) leads to a marginally smaller value of m^, and the choice mq/mq = 120 (lower 
right) leads to significantly lower values of for any fixed value of The effect of 

including bound-state effects is to increase the range of compatible with the measured 
value of ~ “^q/^g = l-lj decreasing to ~ 20% for mq/mq = 10 to 50. 

Finally, we show in Fig. [^the value of at the endpoint of the coannihilation strip when 
Am = 0 and Q^h'^ = 0.1193 ± 0.0042 (green band), as a function of mq/mq: the brown and 
red contours are for Q^h"^ = 0.05 and 0.15, respectively. The band and contours exhibit the 
inverse of the behaviour of the relic density seen previously in Fig. The neutralino mass 
at low mq/mq is below the maximum value of m^, which has a shallow maximum around 
mq/mq = 10 to 50, and falls sharply when mq/mq > 100, reflecting the effect of a breakdown 
in ^ — X conversion. We conclude that, within the framework studied here, m^ < 8 TeV 
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(rising to ~ 9 TeV when allowing for a factor 2 uncertainty in the bound-state formation 
rate) in the Bino LSP case. 



Figure 6: The value ofm^ at the endpoint of the gluino coannihilation strip when Am = m^ — 
= 0 in the Bino LSP case, as a function of mg/mg. The drop at small m^/mg is due to 
the cancellations between the s-, t- and u-channel diagrams for gluino pair annihilation, and 
that at large mq/mg is due to the decoupling of the gluino and neutralino densities. The green 
band corresponds to the current 3-a range of the dark matter density: = 0.1193±0.0042, 

and the brown and red contours are for = 0.05 and 0.15, respectively. 


6.2 Wino LSP 

We now consider the case of a Wino LSP. The left panel of Fig. [^displays the gluino-Wino 
coannihilation strips for = 0.1193 ± 0.0042 for mq/mq = 10, using the same colour 
codings as for the Bino case (red with neither the Sommerfeld enhancement nor gluino 
bound states, orange including the QCD Sommerfeld enhancement but again no bound- 
state effects, black with both effects included, and purple with the bound-state formation 
rate enhanced by a factor 2). We see that in this case the black coannihilation strip extends 
to ~ 7 TeV. Note that the curves appear to diverge at low m^. The reason is that even in 
the absence of coannihilation (large mq — m^), Wino-Wino annihilations are strong enough 
to suppress the relic density below the density indicated by Planck and other experiments 
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Figure 7: As in Fig. and[^ but for a Wino LSP. 


when < 3 TeV The right panel of Fig. [^shows how at the endpoints of the strips 
varies with m^. As previously, the colours of the lines correspond to the colours of the strips 
in the left panel of Fig. We see that the black line crosses the horizontal green band where 
= 0.1193 ± 0.0042 for ~ 7 TeV. 

The left panel of Fig. is the analogue of Fig. |^for the case of a Wino LSP, with the 
green band corresponding to = 0.1193 ± 0.0042 and the brown and red contours to 
= 0.05 and 0.15. We see that is within the preferred range for ~ 7 TeV over 
a broad range 5 < mq/nig < 100. The percentage increase in the allowed range of due to 
bound-state effects, as a function of mg/mg, is similar to the Bino case. As shown in Fig. 
for the Bino case, the fall in the to lower values of is due to the breakdown of g — x 
conversion. The curve hits a plateau for mg/mg > 300 which represents the decoupling limit 
at ~ 3 TeV. 


6.3 Higgsino LSP 


We now consider the case of a Higgsino LSP. The left panel of Fig. displays the gluino- 
Higgsino coannihilation strips for = 0.1193 ± 0.0042 and mg/mg = 10 using the same 


^The Sommerfeld enhancements in the calculations of the thermal relic abundance of a Wino- or Higgsino- 
like LSP were discussed in detail in 34 . We have included an estimate of the Sommerfeld enhancement 
factor for the Wino-Wino annihilations, and we get a similar result as 16 , namely ~ 3 TeV, in the limit 
where the effect of gluino coannihilation is absent. We do not include the Sommerfeld enhancement for the 
Higgsino-Higgsino annihilations, since its effect is much milder in the relic abundance calculations, as shown 
in the literature. 
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Figure 8: As in Fig. but for a Wino LSP (left panel) and a Higgsino LSP (right panel). 


colour codings as for the Bino and Wino cases (red with neither the Sommerfeld enhancement 
nor gluino bound states, orange including the QCD Sommerfeld enhancement but again 
no bound-state effects, black with both effects included, and purple with the bound-state 
formation rate enhanced by a factor 2). In this case the black strip extends to ~ 6 TeV at 
the endpoint where Am = 0. The Higgsino couplings depend on tan (3 and in our calculations 
we have taken tan (3 = 10. Our results are very weakly dependent on this choice. Once again 
we see a divergence of the contours at low m^. In this case, when m^ < 1.2 TeV, Higgsino- 
Hiiggsino annihilations are sufficient to reduce the relic density below the Planck density. 
The right panel of Fig. |^shows how at the endpoints varies with with the colours of 
the lines corresponding again to the colours of the strips in the left panel of Fig. The black 
line crosses the horizontal green band where = 0.1193 ± 0.0042 for m^ ~ 6 TeV. As 

seen in the right panel of Fig. similar values of are found for a range 5 < mq/nfig ^ 100, 
with the drops in the contours to lower values of again being due to cross section 
cancellations at low mq/irtg and due to the breakdown oi g — x conversion at high mq/rug. 
As in the case of the Wino, the curves drop to a plateau for mq/rUg > 300 representing the 
decoupling limit. In this case, the asymptotic value of m^ is ~ 1.2 TeV. The percentage 
increase in the allowed range of due to bound-state effects is again similar to the Bino 
case. 

The decreases in the maximum values of m^ allowed in the Wino and Higgsino cases. 


compared to the Bino case, are due to the effect noted in 35 , namely that coannihilations 
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Figure 9: As in Fig. and[^ hut for a Higgsino LSP. 

may, under some circumstances, increase the relic abundance by coupling ‘parasitic’ degrees 
of freedom. In the Bino (Wino) (Higgsino) case, there are 2 (6) (8) electroweak degrees of 
freedom, linked by coannihilation to the gluinos, that contribute incrementally to the relic 
abundance. This effect is compensated by the decreases in the maximum values of that 
we hnd in the Wino and Higgsino cases. 

7 Summary 

We have studied in this paper MSSM scenarios in which the LSP is (almost) degenerate 
with the gluino, exploring the characteristics and locating the endpoints of the gluino-LSP 
coannihilation strip in the cases where the LSP is the Bino, the neutral Wino or a neutral 
Higgsino. Important ingredients in our analysis are the Sommerfeld enhancement of gluino 
annihilation rates, gluino-gluino bound-state formation and gluino-neutralino conversion. As 
we show, these can affect signihcantly the preferred range of the gluino-LSP mass difference 
along the coannihilation strip, and also the position of the endpoint. 

In the Bino LSP case, we hnd that at the endpoint the LSP mass ~ 8 TeV, increasing 
to ~ 9 TeV if we allow for a factor 2 increase in the bound-state formation rate above our 
calculations. These values are decreased by ~ 1 TeV if the LSP is a Wino, and by a further 
~ 1 TeV if it is a neutral Higgsino. The upper limit on the LSP mass of ~ 8 TeV is weakly 
sensitive to the squark mass for 10 < mq/rug < 50, but is substantially reduced for either 
smaller or larger values of mq/nig. In all cases, the percentage increase in the allowed range 
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of due to bound-state effects may be as large as 50%. 

We are loath to claim that our upper limit on the LSP is absolute, but we do note that 
it is substantially higher than what is possible along the stop coannihilation strip, reflecting 
the larger annihilation rates that are possible for the gluino because of its larger colour 
charge. However, these annihilation rates also depend on the masses of other sparticles, 
notably the squarks in the gluino NLSP case studied here. As have shown, the decrease in 
upper limit on the LSP mass for small mq/rrig is due to cancellations in the annihilation 
matrix elements, whilst the decrease at large mq/rrig is due to the breakdown of glnino- 
LSP conversion. However, we have not stndied the limit mq/rrig —)■ 1, where many more 
coannihilation processes wonld come into play, as might also be the case in non-minimal 
snpersymmetric models. 

Nevertheless, onr analysis does show that a large mass reach to at least 8 TeV will be 
necessary to explore conclusively the possibility of supersymmetric dark matter within the 
MSSM and a conventional cosmological framework. 
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Appendix A 


In this appendix, we hrst recall our procedure for computing the thermal-averaged velocity- 
weighted cross-sections {ai 2 Vre\) for the process 1-1-2—)-3-|-4asis necessary for solving the 
Boltzmann eqnations in Section 5 in an efficient manner, extending the approach used in 
More details of our approach can be found in 
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We start with the sqnared transition matrix elements |Tp (snmmed over hnal spins and 
colors, averaged over initial spins and colors) for the coannihilation processes of interest, 
which are given here in Appendix B, expressed as fnnctions of the Mandelstam variables s, 
t, u. We then express |Tp in terms of s and the scattering angle in the centre-of-mass 
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frame, as described in 36 . Next we define 

' '' 


wis = - 


(fp4 


+P2-P3- Pa) | 7 "p 


327r sV2 


(27r)32E3 (27r)32E4 

P3{s)_ 


dcos^cM |Tp 


( 42 ) 


The total annihilation cross section (Ti 2 (s) is given in terms of w{s) by (Ti 2 (s) = w{s)/s^^‘^pi{s 
The usual partial-wave expansion can be obtained by expanding |Tp in powers of pi(s)/mi 
The odd powers vanish upon integration over ^cm, while the zeroth- and second-order terms 
correspond to the usual s and p waves, respectively. We can therefore evaluate the s- and 
p-wave contributions to w{s) simply by evaluating |Tp at two different values of cos 6 *cm- 


The proper procedure for thermal averaging has been discussed in 37,41 for the case 


mi = m 2 , and in 38,40 for the case mi 7 ^ m 2 , so we do not discuss it in detail here. One 


finds the coefficients a and b in the expansion ([^ of the thermal-averaged cross-sections for 
the processes of interest: 


(ffi 2 t'rei> = an + l>nX + 0(x •') , 


(43) 


where x = rrii/T (assuming mi < m 2 ) by following the prescription given in 36 , using the 
transition amplitudes listed in Appendix B for each final state. When the conversion rates 
are large compared to the Hubble rate, these amplitudes can be used to compute the total 


effective coefficients Ues and b^s by performing the sum over initial states as in (38), and we 


then integrate the rate equation (35) numerically to obtain the relic density. 


Appendix B 


The 2 —)• 2 (co)annihilation processes relevant to the gluino-neutralino (and, in the Wino 
and Higgsino cases, gluino-charginos) system are gg —)■ gg, gg —)■ qAQB, gXi qAqB and 
gx^ —t UAds, where the indices A,B = 1,2,3 for three generations, the neutralino index 
i = 1, ...,4 and the chargino index j = 1, 2. 

Here we list the |Tp for each of these processes, separating the contributions from s-, 
t- and M-channel diagrams. In the following expressions, hnal spins and colors are summed 
over, and initial spins are averaged over. A factor Cini is used to average over initial colors. 
Therefore, |Tp takes the form 


1 — Cini{Ts'^'Ts -\- Tt^Tt Tu^Tu Ts^Tt -\- Ts^Tu Tt^Tu) ■ 


^Our w{s) is also the same as w{s) in 


37 


39 


which is written as W/A in 


40 


( 44 ) 
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We note that there is also the charge-conjugated process for the chargino, gXj —t UAds, 
which we do not list separately. 


gg gg 


There is an s-channel gluon-exchange diagram, and t- and w-channel gluino-exchange dia¬ 
grams. We note that, because there are two identical gluons in the hnal state, an extra 


factor of 1/2 is needed when performing the momentum integration in (42). 
We hnd 


Qm — 
'Ts^'Ts — 
TtxTt = 


1 

M ’ 

11527r^a^ 


-{t- uf] , 


23047r2a2 


'TsX'Tt — 
TtxTu = 


{ml-ty 


57671^0^ 


m? [s‘^{t + 3u) + 2s(t^ -I- 2u^) + 2{t + 

-I- — 2s{t -I- 2u) — 6{t + -|- 2m| [s + 4{t + u)] 

— 4m| — tu -I- 2su + 2{t^ -|- u^)\ | , 


r\ \s{t - u) (4m? -t + u) d- + s^{u -t) + {t - 
(m| — tj 


23047r2a2 


[mt — t) (m? — m) 


(m| - tu) [-4(f -1- M)m? -F 8m| + {t - uf) , 


9 J \ 9 

and TgXTu and TuXTu are related to 7^ x 7/ and TtXTt, respectively, by exchanging t u 
in the corresponding expressions. 


gg QAqB, gx^ gx'j UAds 

These three processes all have t- and w-channel squark-exchange diagrams, and gg -A Qa^b 
also has an s-channel gluon-exchange diagram, whereas the other two processes do not (hence 
TsxTs = TsxTt = Ts xTu = 0 for them). Apart from the couplings, these three processes have 
the same structures as for 7/x7/, TuXTu and TtxTu- In the case of gXj -A UAds) because the 
quark CKM matrix is involved in the chargino-quark-squark vertex, the indices A and B can 
be different even if we restrict to the case of no generation mixing with only left-right mixing 
in the third generation for the up-type and down-type squarks. Therefore, it is convenient 
to write a 6 x 6 up-type squark mixing matrix, , which relates the interaction eigenstates 
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and mass eigenstates of the up-type squarks as follows: 


/ ul \ /lOOO 0 0 \/«i\ 

Cl 0 0 1 0 0 0 U2 

II _ 0 0 0 0 cos 9i — sin 6i Ci 

Ur 0 1 0 0 0 0 C2 ’ 

cr 0 0 0 1 0 0 ii 

\ / \ 0 0 0 0 sin 6 ^^ cos 9^ / \ h / 


( 45 ) 


where 9^ is the stop left-right mixing angle. The mass eigenvalues are correspondingly dehned 
as rufy^ = ma^, = RRcii = RRc 2 -i R^jj^ = R^h mp^ = mi^. A similar 

mixing matrix, , is introduced for the down-type squarks, by changing 9i to the sbottom 
left-right mixing angle, 9^. The mass eigenvalues RRl,^_^ are also defined similarly. 

For gg —)■ Qa^b we hnd 


'TsX'Ts 

'Ts^'Tt 


3847r^(5ABtt= 


\^m\ [2ml + s — t — u) + 2ml(s — t — u) + 2m\ + 2ml + t‘^ + u^ , 


p 1927r^(5ABas 
^ s (t — ml ^ 


mj {2ml + s-2t)+ml{s- 2t) +m\ + ml + t^] + \Z{p 


- 2 mim 3 {3ml + 3^^ -2t-u) (^zQ * + z{^ (^( 4 + 3 )p) * 


where mi = my and m 3 = m/^. The index f = U,D, the index f = U,D. mjy^ 2 3 = mu,c,t, 
mDi 2 3 = RRd,s,bi and 7^ x TJ, is related to 7^ x 7t by exchanging t -H- u. 
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For all three processes, TtxTt, TuXTu and TtxTu take the following forms: 


r.xT, = -W7r-TTP’">“h/i(Ap)A(A9)- + /H(Ap)/i(A<()*) 

P,q=l V ^tp) V 

+ {ml + ml-t) {fi{A,p)fl{A,qy + fji{A,p)fji{A,qy)] 

X [2m2m4 p)gR{B, q)* + g^^{B, p)g\{B, q)*) 

+ {ml + ml-t) [gL{B,p)g\{B,qy + g^p^{B,p)g^p^{B,qy)] , 

T.yT^ = It- 

P,g=i [t-mfj 

X I [fiiA qTfLiA p) {m 2 gUB, g)* [m^ {ml + ml-t) gl{B, p) +mi{ml + ml-t-u) gl^^B, 
+ gL{B, q)* [{m\ml + mlml - tu) gl{B,p) + mim4 {mj + ml - u) gR{B,p)] | 

+ msfliA, qyf^iA, p) |m2^i(5, q)* [{ml + ml-u) gl{B, p) + 2mim^g^AB, p)] 

+ gR{B, q)* [mi {ml + ml-t- u) gl{B,p) + mi {ml + ml-t) g^iB^p)] 

+ (L <H- i?) I + (t <H- M, mi <H- m 2 ) , 

where the (L <H- i?) in 7t x 7)^ applies to all the L and R in the indices, and the {t <H- u) 
applies to both the t and u in the indices and the Mandelstam variables. Again, T^xTu is 
related to TtxTt by exchanging mi <H- m 2 and f <H- m in both the indices and the Mandelstam 
variables. 

The couplings and masses involved in the above expressions are listed below. 



For gg -)■ qAQB- 




c 


tu 


mi 


m3 

IUAp) 

9iiB,p) 

9r{B,p) 


1 


M ’ 

^UU _ 

32 

’ 

256 

IT’ 

m2 = 

m-g, 

7714 = 

rrifA , 

= 


mAp) 

— 7^ 

- ^(A+3)p ’ 

fniAp) 

= -zip. 

9l{B,p) 

= K);- 

9l{B,p) 

= “ (^(S+3)p 


where the index f = U,D, and the index f = U,D. 


29 



For gx° qAQB- 


^ini 

mi 

m 2 

mf 

Lp 

/i(Ap) 

/MAp) 


1 

8 ’ 


^UU _ _^tu _ 


mg, 


vrac 


m^o , 




yf 

^{A+3)p : 

— 7^ 

^Ap > 


W2mf^c\ ( ^/b+3)p 


V2'. 




9l{B,p) 

g^iB^p) = iV2g2NiiQfgtan{e^) [Z^ 

mAp) 


- tV 2 g 2 [ZQ * [(iVn) * tan (6^ {Qf, - T|J + T|^ (iVi^) *] , 

A2mf^A i^ip] 


{B+3)p 


V2: 


m,. 


• /o n (AT \ * I (a \ yf '^92fT^fA^L^Ap 

iA/2g2QfA{Nii) tan{e^)ZI 


{A+3)p 


a/2\ 


fniAp) 

9l{B,p) 

9l{B,p) 


»92’»/AZ(.1+3). _ [JV.J J3J _ 

V ^’((T'w 


yf 

^Bp 


^ ^{B+3)p 


where N is the 4x4 neutralino mixing matrix as dehned in 42 , and g 2 is the Standard 
Model SU(2)^ conpling constant. For np-type qnark hnal states, the index f = U, and 

1 


rTi6 

^Ia 

= Tl = 

2 ’ 

QfA 

= Qjb = 

2 

3 ’ 

4 

= A = 

(4) 

m 3 

= rnuA, 


m4 

= fnug 3 


^fA 

= mUA, 


^fB 

= rnug . 
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taii/S = V 2 /V 1 , and the vacuum expectation values of the two Higgs doublets are dehned as, 


For down-type quark dual states, the index f = D, and 

rp3 _ rp3 _ ^ 

^fA - -^fB - 2 ’ 

Qia ~ Qib ~ ~3 ’ 

4 = = (4)* = (4)- = sec(,3) (/V„) •, 

ms = mos, 

^4 = m£,g , 

^fA = 

TTIJb ^Dg ■ 
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For gx\ UAds- 


Cin.i ^ 


1 

8 ’ 


c** 

mi 

m 2 

m 3 

7714 

7nu„ 


^UU _ 


-c*" = 


TT 0.0 


= m~g, 

= 

'n^Dg ) 




= rrif 


fiiAp) 

fniAp) 


yU 

^(A+3)p 5 

— 7^ 

^Ap y 


9iiB,p) = 


ig2 csc{/3)KcBmuc (^ 2 ) ^ 


^(C+3)p 


9UB,p) 

mAp) 


C=1 
3 

E 

C=1 
3 

C=1 \ 


\/27nu 

ig2 sec{/3)KcBUj2mDi: 


— i92KcB Ai] 


yU 

^Cp 


yU 

^Cp 


A/2m^ 

I ig 2 KAc csc(/ 3 )mc/^ (1^2) 


¥ rvD 
^Cp 


72. 






V \/2m^ ‘ y 

^7L(i^,P) = 

where the K matrix is the quark CKM matrix, and U and V are the 2x2 chargino mixing 


zi 




matrices as dehned in 42 


Finally, we give the s-wave result (i.e., the coefficient a in Eq. (43)) for the gg —)• qAQB 
channel, in the limit of a common squark mass and massless quarks, with no generation or 
left-right mixing in the squark mixing matrices (the case considered in the main body of the 
text). In this limit, the contributions from each of the six quark flavor final states are the 
same, and the result of putting all the six quark flavors together is 


^gg^QA^B limit value — 


dual ( 


rrig 


m 


i) 


8m? (m? -I- m?) ^ 


(46) 
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where rriq is the common squark mass. When rriq S> rrig, only the s-channel gluon-exchange 
diagram contributes, and the above expression is proportional to On the other hand, 

when rriq —)■ rrig, the above expression approaches zero. This cancellation of the s-, t- and 
M-channel contributions results in the feature of the plots at small values of rriq/rrig that are 
commented upon in the main body of the text. 
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